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Abstract. A probabilistic cellular automaton for cargo transport is pre- 
sented that generalizes the totally asymmetric exclusion process with a 
defect from continuous time to parallel dynamics. It appears as an under- 
lying principle in cellular automata for traffic flow with non-local jumps 
for the kinetic constraint to drive as fast as possible. The exactly solvable 
model shows a discontinuous phase transition between two regions with 
different cargo velocities. 
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1 Introduction 

Non-equilibrium phase transitions can rarely be calculated exactly, i.e. without 
need of approximations or fits of numerical data. One paradigmatic system where 
this is possible is the totally asymmetric simple exclusion process (TASEP) (see 
[T] and references therein). The model is defined on a Id discrete lattice with 
sites being either empty (holes) or occupied by a single particle. A randomly 
chosen particle moves to the right at rate 1 provided that the site is empty. For 
open boundaries with particle input at rate a and output at rate /3 this leads 
to three different phases: a low-density, a high-density and a maximum-current 
phase. For finite systems the process can be solved exactly by the matrix-product 
Ansatz (see [2] for a recent and exhaustive review) and the complete thermo- 
dynamic behavior that is relevant for understanding the phase diagram can be 
extracted from the asymptotics of this solution. On the ring the process has a 
uniform groundstate pQ. However the presence of a defect particle leads to a 
rich phase behavior [3j. In the defect TASEP usual particles move 10 — > 01 at 
rate 1, the single defect particle moves itself forward 20 — > 02 at rate a and 
can be overtaken by usual particles 12 — > 21 with f3 The solution is formally 
related to the open-boundary case. There is one shock phase and three phases 
where it behaves once like a particle, once like a hole and once like a second-class 
particle. The second-class particle case corresponds to a = /3 = 1. To the left 
the second-class particle looks like a hole (as seen from a particle) and to the 
right it looks like a particle (as seen from a hole). This case was studied in @] 
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since it can be used to study the shock in the TASEP on the infinite line with 
step- initial condition. The connection is that, due to its very special dynamics, 
the defect 2 samples the random shock position on preferring configurations like 
000021111. 

The defect TASEP can alternatively be understood as a cargo transport process: 
The defect is a usual particle carrying a cargo which can be handed over to a 
particle behind. See [5] for a different definition of cargo in the same context. 
These mechanisms play a fundamental role in the biology of intracellular trans- 
port. Processive motors like kinesin transport cargo over long distances [6]. The 
cargo can alternatively be interpreted as virus particles that use carrier particles 
in order to attain the interior of a cell, see [Fj for further references. 
The Nagel-Schreckenberg model [7] is often referred to as the minimal model for 
one-lane traffic-flow on a freeway. There it is essential to allow for faster and 
slower cars to get a realistic flow-density relation: cars can move up to v max 
sites per time step. Further all cars are updated simultaneously according to a 
parallel update and move independently with probability p. For w max = 1 it is 
equivalent to the TASEP with parallel dynamics. The steady state on the ring 
shows nearest-neighbor correlations and has a simple pair-factorized form [8]. 
For open boundaries the matrix-product technique could be generalized to ob- 
tain the exact steady state |9ll0j . It can be interpreted as a pair-factorized state 
as on the ring modulated by a matrix-product state |llj . 

In this article we introduce a generalization of the cargo-transport process to 
discrete time with parallel updating and give its exact solution. Here we restrict 
ourselves to light cargo, i.e. the case where the speed of a particle is not low- 
ered by the presence of cargo. This appears naturally in the steady state of a 
traffic cellular automaton [TT]. We will see that non-local jump processes where 
particles drive as fast as possible can lead for non-deterministic hopping to a 
discontinuous phase transition on the ring. 

2 The Cellular- Automaton Model and its Solution 

Consider a periodic one-dimensional lattice with sites being either occupied by a 
particle (in state r = 1) or empty (r = 0). One of the particles carries a cargo to 
which we refer to as a defect (r = 2). The particles are updated simultaneously 
and every particle (with or without cargo) moves forward with probability p. If 
the site behind it is occupied, the cargo carrying particle can independently give 
its cargo back at probability /3. This simple dynamics is encoded in detail in the 
transitions 



10 — ^ 01, at rate p, 
020 -> x02, at rate p, 
120 -> 210, at rate (3(1 



—> 201, at rate /3p, 
121 -> 21x, at rate /3, 



-> 102, at rate (1 



P)p, 



(1) 
(2) 
(3) 
(4) 
(5) 
(6) 
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with x being either or 1 indicating that the site can be either empty or occupied 
due to the parallel update. For example the evolution of the pattern 020 can be 
affected by a particle to the left moving itself forward. This is quite a general 
scenario that might apply to biological intracellular cargo transport. The parallel 
update reflects highly active transport where many particles move at the same 
time. 

For j3 — the cargo is attached to one special particle for all times and its 
dynamics is the same as for the other particles. This corresponds to the usual 
TASEP (with parallel update) and the single occupation r = 2 can be replaced 
by r = 1. The steady state for a lattice with L + 1 sites (with one of them 
occupied by the defect) is 

L+l 

P( Tl ,T 2 ,...,T L+1 )= Y[P({Ti-i,Ti}), (7) 
i=l 

thus it factorizes [5] into symmetric two-site factors P{ri-\Ti) = P({r,_i,rj}) 
with 

P(00) = 1 - p - J/p, (8) 
P(10) = J/p, (9) 
P(ll) = p-J/ P . (10) 

Here J is the particle current 

We found [IT] that the steady state for (3 > can be calculated exactly too. 
Here §7§ is generalized to 

P(2,T 1 , ...,T L ) (X /(TO/faTs) • ..f{T L -lT L )f{T L ) 



(W\ 



Y[ Ti D + (1 - Ti )E 



i>l 



\V) (12) 



This is a pair-factorized state (mainly the steady state for j3 = 0) modulated 
by a matrix-product state. Up to the normalization this is very related to the 
TASEP with open boundaries |9lllj . The vectors (W\ and \V) represent the 
defect and the matrices D and E represent particles and holes respectively. The 
operators obey the algebra 

(W\EE = (1 -p)(W\E, (13) 
(W\ED = (l-p)({W\D+p), (14) 
DE = (1 — p)[D + E +p] , (15) 

D\V) = ^j^\V), (16) 
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for details see [II]. Note that the relations (15) and (16) are quadratic and the 
other relations are cubic. Accordingly the dynamical rules in (T|6) are quadratic 
(for 10 and 12) and cubic (for 20) respectively. The thermodynamic particle 
current (11) obviously is unaffected by the cargo. 



This steady state appears also in cellular-automaton models for traffic flow with 
non-local jumps under kinetic constraint. Consider the following process on a 
periodic one-dimensional lattice with sites being either occupied by one car (in 
state r = 1) or empty (r = 0). The update rules applied simultaneously to all 
cars (= particles) are 

100 -> 001, with probability p, 

101 — > Olx, with probability /?, 

where x denotes either a particle or hole. The maximum velocity i» max thus is two 
sites per time step instead of one in the usual TASEP and the kinetic constraint is 
that cars can not drive at reduced speed 1 if they could move at maximum speed. 
This leads under the parallel dynamics to a non-local repulsion between cars, so 
that finally in the thermodynamic limit only even gaps (0,2,4, . . . holes) have 
non- vanishing probability. Figure [I] a) shows schematically the allowed moves in 
a stationary configuration. For even number of holes the process is equivalent to 



a) f' ? 




•ob 



o< 



OO 



Ti Ti+l 



Fig. 1. a) Allowed moves in a stationary configuration: only one odd gap between 
particles b) Equivalence of moves and according reduction of lattice units 



the TASEP and for odd number of holes it is equivalent to the cargo-transport 
process. In this case a single hole in an environment of particles and hole pairs 
is formed that plays the role of cargo attached to varying particles. In figure [T] 
b) one sees that the particle movement of a single site is equivalent to backward 
movement of the 01 position. The 01 pair plays the role of the defect, having the 
characteristic 'Janus face', looking to the left like a hole and to the right like a 
particle. Usual holes are replaced by hole twins 00. To be precise, the probabilities 
( 8p| would be rewritten here as P(00) = P(0000), P(01) = P(001) = P(0001). 
In the following we restrict ourselves in the terminology to the cargo-transport 
process. 
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3 The Phase Diagram 

There is a discontinuous phase transition at the density 

0(1-0) 



(17) 



For p < p c and p > p c one hnds different velocities of the defect which can be 



Q.O. 




Fig. 2. The phase diagram shows the p-(3 plane for p = 3/4. The thick line 
separates the two phases and on the dashed line the velocity of the defect changes 
its sign. 



calculated through 



v=p(l-p + )(l-/3p-)-l3 P -. 



(18) 



Here p- and p+ are the densities directly behind and in front of the defect 2. 
The dynamics of the defect is obtained from ( l|6 1 . It moves either forward with 
probability p if it has a hole in front while at the same time there is no particle 



directly behind that simultaneously catches the cargo: first term in (18 1. Or it 



moves backwards with probability /3 if it has a particle behind: second term in 



( 18 ). The neighboring densities are in terms of the current J(p) defined in ( 11 ): 



1 



J\p-Jf 



P- 



P{\-pYp{ p - j) 2 + (i-p)j 2 

PR - J 



. PR ~ PJ 
J 

PR, 

P-P J 
,p 2 (l-(3)p 



, for p < p c 
for p > p c 



for p < p c 
for p > p c 



(19) 



(20) 
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Note that there are many ways to express the results due to the relation J(l — 
J) = pp(l — p). In other words the square root in J appears in every power of J 
with a certain prefactor. Equations ( 19|20 ) yield 



P 



l-2p 

1 - 2J' 
J-Pp 

.vp- PJ 



for p < p c 
for p> p c 



(21) 



The defect velocity vanishes for 

P~P 



Po 



P~P 2 
1/2, 



for P < 1 - \/l -p, 



(22) 



for P > 1 - y/1 -p, 



and is positive for p > po and negative for p < p$. This leads to the phase 
diagram, depicted in figure [2j The formulae §Tf\ and ([22]) can alternatively be 




P 

Fig. 3. Velocity of the defect for p = 3/4 for varying /3. Continuous line: result 
for p < p c (P)- The broken lines correspond to P = 0, P = 0.25, /? = 0.5, /3 = 0.67 
from top (dashed) to bottom (dotted) and are valid for p > p c (P)- The segmented 
line is the velocity of holes 



interpreted in terms of a critical value p c = P(p c ) and Po = P(po). This gives 

p c = J(p c )/(l- Pc ), Po = J(p c )/ Pc . (23) 

The value of p c and /?o respectively then is essentially the absolute velocity of 
holes and particles at the transition. Figure [3] shows the character of the defect 
and the discontinuous phase transition. For P = its velocity is given by the 
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upper curve and equals the velocity of particles. In the zero-density limit the 
velocity is independently of ft equal to p. For increasing ft the second phase 
appears: The velocity of the defect jumps at p = p c (ft) from the continuous 
curve to the corresponding dashed curve. p c increases with ft until ft = p where 
p c = 1, so that the system is completely in the second phase for all densities. 
Note that for ft = 1/2 one has p c = 1/2. Finally in the limit of the fully occupied 
lattice p = 1 and p c < 1 the velocity equals ft for ft < p. However for p c > 1 
it can never increase the value of p which there is the corresponding velocity of 
the holes given by the segmented line. 

For the phase p < p c , which is purely present for ft > p, it is important to stress 
that the quantities p+, ftp- and v are independent of ft. The density profile is 
symmetric around the defect and has an algebraic decay. This corresponds to the 
second-class particle phase in the defect TASEP mentioned in the introduction. 
As in continuous time the defect velocity (21 1 is given by v = dJ /dp which has 
the form of a group velocity and becomes v = 1 — 2p for small p, compare table 
[I] Thus the defect travels with the velocity of the density disturbance. 



4 Limits 

Table |l]) shows the limit of small hopping probabilities: p = dt, ft = ftdt with 
dt — > 0. Note in comparison that the velocity of normal particles is always J / p 
which gives p(l — p) +p 2 (l — p) 2 + . . . This limit has been studied in the traffic 
picture in |12ll3j and corresponds to the defect TASEP. As mentioned before, 



p 


V 


l-p+ 


P- 


P< pc 
P> Pc 


(1 - 2p)dt 
p)dt 


(1 - pf [1 + 2p(l - p)dt] 
(l-P)(l-p) [l + + p(l-p))dt] 


p 2 /ft [1 - (1 - P )(l - 3p)dt] 
p[l-(l-p)(l-ft-p)dt] 



Table 1. Continuous-time limit p = dt, ft = ftdt: Values to the order 0(dt) of 
velocity, empty-space density in front and particle density behind the defect in 
the two phases. Here p c ~ ft [1 — (3(1 — ft)dt] 



in the limit (3 = the defect moves only forward and loses its role as a defect. 
Therefore the steady state is the same as for the TASEP. Thus one has the 
expressions given in table (pj. In comparison the results from (|19|20|21| for 



7WV 



p+ 



J_ 
pp 



i-L 
pp 



Table 2. Limit of the discrete-time TASEP ft = 
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p > p c are rewritten and expanded around the TASEP value: 



P-(P) = 
l-p+{P) = 
v((3) = 



PUD 



P(ll) 



l 



P-PP(10) P 
p-P P(10) _ P(10) 

J-/?p 



P(10) 



p 



p-pp(io) p V pp 2 



(24) 
(25) 
(26) 



One sees that p_ is mainly the same as for /3 = but the density to which the 
numerator is adressed is reduced by backward moving so that /?_ is increased. 
The same holds for the velocity v. 1 — p + is even the same as for /3 = up to a 



scale which is, using (17) and (22), given by p c /po- 

For p = 1 (24)-(26) yield the results displayed in table[3j For p < 1/2 all particles 



p 


V 


1-P+ 


P- 


p< 1/2 
p> 1/2 


1 - 


1 

"P- 


/3p 


1 

1-P 




2p- 1 


P~ 


0(1 


-P) 


P 


P 


-P{l-p) 



Table 3. The partially deterministic case p = 1: Velocity, empty-space density 
in front and particle density behind the defect. 



are separated and move deterministically as in the TASEP. For p > 1/2 the effect 
of P on the velocity remains for all possible values and the phase transition 
disappears. 



5 Conclusions 



A cellular automaton for cargo transport was introduced that generalizes the 
(continuous-time) defect TASEP. The parallel update is often more realistic in 
describing active many-particle transport and makes the link between determin- 
istic and random-sequential dynamics. The point of interest was a single defect, 
i.e. a particle carrying light cargo in an environment of particles and holes on a 
periodic Id lattice. Particles move forward with probability p and if a particle is 
directly behind the particle that carries the cargo it may catch the cargo with 
probability p. We found a discontinuous phase transition between two phases 
with different cargo velocities. Successively increasing p lowers its velocity only 
until p — p. Then a saturation effect appears where the velocity becomes inde- 
pendent of p. The same holds for the stationary state of a cellular automaton 
for traffic where the 'cargo' corresponds to small headway that is formed dy- 
namically. It is attached to one car from behind. If the subsequent car comes 
close enough it will catch it up. The fact that the cargo process appears in a 
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seemingly unrelated non-local jump process underlines its universal role. The 
exact matrix-product state and its cubic algebra holds also for the two-species 
case where multiple cargo is present. It is also trivially generalized to case where 
cargo lowers the speed of the particle to a < p. For p — j3 and in the presence 
of several second-class particles it serves also as a model system for the study of 
shocks on the infinite line. 
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